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We study the carbon-dope aluminum clusters by using 
time-of-flight mass spectrum experiments and ab initio cal- 
culations. Mass abundance distributions are obtained for an- 
ionic aluminum and aluminum-carbon mixed clusters. Be- 
sides the well-known magic aluminum clusters such as Al^ 
and Al^" 3 , AI7C - cluster is found to be particularly stable 
among those Al n C~ clusters. Density functional calculations 
are performed to determine the ground state structures of 
Al n C~ clusters. Our results show that the Al7C~ is a magic 
cluster with extremely high stability, which might serve as 
building block of the cluster-assembled materials. 
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In recent years, clusters and cluster-based materials 
have been a field of intensive research,dj|ie to both funda- 
mental and technological importanceul3. The structural, 
electronic, magnetic, and optical properties of the clus- 
ters are different from those of constitute atoms or bulk 
phase and depend |-seirsitively on the size and composi- 
tion of the clusteruS. It is desirable to assemble the 
cluster-based materials from properly designed clusters 
so that the unique, properties of these individual clusters 
can be retainedoUJ. To be a building block of cluster- 
assembled materials, the cluster should be highly stable 
and relatively unreactive. Thus, the clusters would inter- 
act weakly with each other and maintain their identities 
when they are brought together in the cluster-assembled 
solids. A well-known example is the Cgo solidail Besides 
Ceo, it has also been suggested that the metal clusters 
with electronic and geometric shell structure can be used 
to construct cluster-assembled materialsEl. In this direc- 
tion, some previous theoretigaLpffieirts have been devoted 
to aluminum-based clustersQ'QoEj. First principles cal- 
culations predicted that doping in the aluminum clusters 
can enhance the stability of certain magic clusters like 
AI13 and modify the physical and chemical properties of 
the clustersBO. The stability and electronic properties 
of AI13K solid has also been calculated^. 

Experimentally, particular attention has been paid to 
pure and carbon-doped aluminum clusters. The mass 
spectra of pure and doped-AL clusters have been ob- 
tained by different groupsc-atll. L.S.Wang et al. have 
performed a combined photoelectron spectroscopy and ab 
initio study on the small neutral and-anipnic aluminum- 
carbon clusters A1„C~ (n — 3 — 5)c3~El The reaction 
of Al n C~ clusters with oxygen was examined by CastLe.- 
man's group and AI7C - was found as a magic clustered. 



In this work, we shall investigate the structure and stabil- 
ity of A1„C~ clusters via a combination of experimental 
cluster mass spectroscopy and ab initio calculations. 

Our experimental apparatus consists of a standard 
Smalley-type laser vaporization/molecular beam cluster 
source, and a Wiley^McLaren time-of-flight mass spec- 
trometer (TOF MS)E3. We use a Q-switched frequency 
doubled Nd: YAG laser [532nm, 15ns full width at half 
maximum (FWHM)] to vaporize the 6 mm graphite tar- 
get rod. The aluminum target contains a trace amount 
of carbon with C/Al ratio as 3xl0 -4 . The laser spot on 
the target is less than 1 mm and the typical laser output 
is about 40 mJ. Pulsed He carrier gas stream (purity: 
99.995%) generated from a pulsed valve is used to cool 
down the plasma above the target that has been pro- 
duced during the laser vaporization process. The clusters 
synthesized in the nozzle channel will then be carried into 
vacuum via supersonic expansion. The high pressure of 
the He carrier gas is necessary for generating the anionic 
Al~ and A1„C~ clusters in a mass range from n = 3 
through n = 38. The clusters form a cluster beam af- 
ter passing a skimmer (2 mm, Beam Dynamics) located 
about 2 cm downstream from the nozzle exit. 

The cluster mass abundance distribution is analyzed 
by the Wiley-McLaren TOF MS driven with a negative 
high voltage pulse generator (EG & G Model 1211, 100 
ns rise time, -2 kV voltage). The extraction and accel- 
eration regions of the TOF MS are 2 cm and 1 cm long 
respectively. The length of the free-drift tube in the TOF 
MS is 105 cm. An X-deflector is applied immediately af- 
ter the acceleration region to compensate the transverse 
energy of the ionized clusters. By adjusting the volt- 
age exerted on the X-deflector, clusters in different mass 
ranges can be recorded. The clusters are detected by a 
dual-microchannel plate (DMCP) and the signal from the 
DMCP detector is fed into a LeCroy 9350AL digital oscil- 
loscope (500 MHz bandwidth, 1GHz sampling rate and 
2x2M maximum record length). Typically, each mass 
spectrum is obtained by averaging 200 shots or more in 
our experiments. 

In this work, all electron density functional calcu- 
lations on A1„C~ (n=l-13) cluster anions, have been 
performed by using the DMol programed A dou- 
ble numerical basis including d-polarization function 
(DND) are chosen. The density functional is treated 
by the generalized gradient approximation (GGA) with 
exchange-ccpjelation potential parameterized by Wang 
and PerdewEJ. Self-consistent field calculations are done 
with a convergence criterion of 10~ 6 a.u. on the total 
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energy and electron density. For each cluster size, we 
started from a number of possible configurations and per- 
form full geometry optimizations with Broyden-Fletcher- 
Goldfarb-Shanno (BFGS) algorithm. We use conver- 
gence criterion of 1CP 3 a.u on the gradient and displace- 
ment, and 10 a.u. on the total energy in the geometry 
optimization. For the smaller clusters, different possible 
spin multiplicities are also tried. 
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FIG. 1. Mass spectrum for 

aluminum and aluminum-carbon mixed cluster anions, Al~ 
Al n C~ and A1„C^", in the small size range (n = 5 — 17). The 
spectrum peaks for A1„C _ anions and magic Al^" 3 clusters 
are marked. Among the Al n C _ clusters, a global maximum 
at AI7C - and a less notable local maximum at AlgC - can be 
identified. The detailed spectrum structures around AI7C - 
are further illustrated in Fig. 2. 
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FIG. 2. An expanded spectrum exhibits the more detailed 
structure of Al~, Al n C~ and Al n C^ in the mass range from 
n = 5 to n = 8. 

The spectrum of cluster anions Al~ and A1„C~ in 
Fig.l covers the small cluster size ranging from n = 5 
to n = 17. The spectrum is recorded when a helium 
carrier gas stagnation pressure of 8 atm and a deflec- 
tion voltage of 90 V exerted on the X-deflector are used. 
One prominent feature in the spectrum is the well-known 
magic cluster Aljg, which corresponds to a complete 
icosahedral geometrical shell structured and a 2p elec- 
tronic shell closiog with the total number of valence elec- 
trons n e = 40trffiLj. In Fig.l, three sequential maximum 
peaks centered at AI7C - are also found. For clarity, an 



expanded spectrum in the mass region around AljC~ 
with a better resolution is presented in Fig. 2. AlgC^, 
AI7C - and AI7C.7 are identified as local maxima. It 
is noted that under the above experimental conditions, 
these species are always found as local maxima. There- 
fore, one can conclude from the mass spectrum that the 
A^C - has the highest abundance among the anionic 
Al n C _ clusters. The present results on AUG-T agree well 
with previous experimental mass spectratSEd. 

To investigate the aluminum cluster anions in a larger 
mass range, we use a 10 atm He gas stagnation pres- 
sure and a 270 V X-deflection voltage. The spectrum 
of Al~ and A1„C~ anions recorded in such conditions is 
presented in Fig. 3. In Fig. 3, Al^ 3 is obtained as a local 
maximum. The seventy valence electrons in Al^ 3 coincide 
the 3s.eLec±ronic shell closing under the spherical jellium 



modeocniS. The geometric.-^ 
contribute to its stabilityli3t^ 



icture of Al^ may also 
1 Whether the slightly 
intense ALjg (n e = 106) peak in the spectrum has an 
electronic shell origin remains an issue of discussion. Be- 
sides, a sudden increase in intensity of Al^ 7 and AI17C - 
is found in Fig. 3. Our present mass spectrum of Al~ 
compares well with previous resultsta. However, in this 
study, we would like to pay the particular attention to 
the unusually stable species A1 7 C~. 
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FIG. 3. A mass spectrum of Al^" and Al n C~ in the larger 
size range. 

In this paper, we exploit the lowest energy structures of 
the A1„C~ (n = 1 — 13) by using density functional geo- 
metrical minimization. The obtained ground state struc- 
tures of Al n C~ (n = 1 — 7, 8 — 13) are shown in Fig. 4 
and Fig. 5 respectively. In general, the carbon atom is lo- 
cated in the center of the equilibrium structure of A1„C~ 
clusters. The calculated cluster properties such as atom- 
ization energy and HOMO-LUOM gap are described in 
Table I. Our present results on the cluster energv differ- 
ence AE n agree well with previous calculations^. 

The ground state configurations and bond length pa- 
rameters of the smaller A1„C~ (n = 1—4) compare 
well with previous calculations in Rcf.[13] (Fig. 4). Spin 
triple state is found as ground state of A1C~ and its 
bond length— ja-,1. 89 A, which agree well with previous 
calculationa 13 lrT The minimum energy structure found 
for AkiC - is an isosceles triangle (C2 V ) with bond length 
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1.85 A and apex angle 9 = 110.6°, in good agreement 
with the 1.85 A— and 103 ° obtained from previous ab ini- 
tio calculations^. The lowest energy structure of AbjC - 
is a carbon-centered planar triangular structure with C- 
Al bond length 1 9.1 A , which is close to that of 1.90 A in 
previous studiesE£E3. For AI4C - , the lowest energy con- 
figuration is obtained as a carbon-centered planar trapez- 
iform. The Al-C distances in AI4C - are 1.96 and 2.05 A , 
This structure is consistent with that in Ref. [13] but dif- 
ferent from the symmetric constraint calculation done by 
L.S.Wang et a/.Ej. 
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FIG. 4. Lowest-energy structures for small Al^C - 
(n — 1 — 7) clusters. The atoms with larger radius denote 
carbon atom. 



Table I. Theoretical atomization energy (AE) and 
HOMO-LUMO Gap (A) for anionic A1„C~ clusters with 
n = 1 — 13. AE n are the differences in atomization en- 
ergy between cluster with n and n — 1 atoms for A1„C~ 
clusters, while AE n a are those taken from Ref. [13]. 
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AE (eV) 


AE n (eV) 


AE n a (eV) 


Gap (eV) 


AKT 


4.80 






0.62 


AI2C- 


9.20 


4.40 


4.46 


2.40 


Ai 3 cr 


13.44 


4.24 


3.97 


2.15 


auc- 


15.98 


2.54 


2.62 


1.45 


Aiscr 


18.60 


2.62 


2.38 


0.80 


Ai 6 cr 


21.64 


3.04 


3.21 


0.79 


AiyC- 


25.66 


4.02 


3.91 


1.72 


Aigcr 


27.52 


1.86 


1.94 


0.36 


A1 9 C- 


30.72 


3.20 




1.48 


AlioC" 


32.95 


2.23 




0.37 


AlnC" 


35.68 


2.73 




0.87 


AliaC" 


39.01 


3.33 




0.21 


AlisC" 


41.46 


2.45 




0.52 



The low-energy structure of A1 5 C~ can be considered 
as being a distorted square pyramid face-capped by one 
Al atom. This structure has been obtained as a local 
minimum m previous studvES. In our calculation, it is en- 
ergetically lower than the equilibrium structure (a com- 
pressed octahedron) found in Ref. [13] by 0.09 eV. The 
low-energy structure for AleC - can be constructed by 
face-capping one more aluminum atom on the configu- 
ration of A^C - in Fig. 4. However, a carbon-centered 
triangular prism was found as ground state structure in 
Ref. [13]. The energy difference between these two iso- 
mers is 0.21 eV from our calculation. In the case of 
AI7C - , we find a central carbon surrounded by a seven- 
atom aluminum cage, in complete agreement with previ- 
ous workEJ. Thus, the particular low reactivity of AI7C" 
against oxygencil can be partially understood by the close 
aluminum cage, which can protected the high reactive 
carbon site from approaching O2 molecules. 
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FIG. 5. Lowest-energy structures for larger Al n C~ 
(n = 8 — 13) clusters. The atoms with larger radius denote 
carbon atom. 

As shown in Fig. 5, the low-energy structures of AlgC - 
are based on face-capping of one Al atom on the. AljC~ 
cage, which agree with previous calculations^. The 
structure for Al g C~ can be viewed as a carbon-centered 
octahedron edge-capped by three atoms. Similar to 
A1 8 C~, the low-energy structures of Ali C~ can be con- 
structed by face-capping of three aluminum atoms on 
the cage-like structure of AlyC~. Starting from n = 11, 
the equilibrium configurations of A1„C~ switch to the 
icosahedron-based packing. AlnC~ is a truncated icosa- 
hedron with slight distortion. AI12C - is a perfect carbon- 
centered icosahedron with C-Al distance as 2.56 A . The 
equilibrium structure of AI13C - is a distorted carbon- 
centered icosahedron capped by one Al atom. 

In Fig. 6, we plot the second differences of cluster 
energies defined by A 2 E(n) = AE(n + 1) + AE(n — 
1) — 2AE'(n), where AE(n) is the atomization energy 
of A1„C~ clusters from DFT calculations. In cluster 
physics, the A2-E(n) is a sensitive quantity that reflect 
the stability of clusters and can be directly compared to 
the experimental relative abundance. In Fig. 6, a global 
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maximum at n = 7 and local maxima at n = 3, 9 are 
obtained. The present theoretical results are consistcejt 
with our experiment and previous works quite wellE-3'cII. 
In addition to the most pronounced peak at AI7C - , 
our experimental mass spectrum in Fig.l also show a 
local maxima at n = 9. On the other hand, Ai3C~ 
was found aa.a local maximum clusters in Castleman's 



experiments 



China. 

* Corresponding author: zhaoj@physics.unc.edu 



pq 
< 





AI 7 C" 


AI3C- 


■ 


■ 


/ \ Al = c 






■ 


\ / \ / " — 

■ 



2 4 6 8 10 12 

Cluster size n 

FIG. 6. Second differences of cluster energies A^Ein) de- 
fined as AE 2 {n) = 2AE(n) - AE(n - I) - AE(n + 1 ) as a 
function of cluster size n for n — 2 — 13. E(n) is the atom- 
ization energy of Al n C~ clusters. Maxima at n —3, 7, 9 are 
found. 

From our experimental and theoretical investigation, 
AI7C - is always obtained as a magic cluster with pro- 
nounced high stability among the anionic Al ra C _ clus- 
ters. In addition to the maximum on the size dependent 
cluster energy difference, AE n or A 2 i?(n), A1?C~ also 
possess a HOMO-LUMO gap of 1.72 eV larger than its 
neighboring clusters (Table I). Such large gap should have 
certain contribution to the low reactivity upon oxygen. 
The particular stability of the AI7C - could be useful in 
construct cluster-based materials. For example, it might 
be possible to assemble A\jC~ with alkali metal ions to 
form stable solid. The further calculations on this direc- 
tion is still under way and the results will be described 
elsewhere. 

In summary, we have performed both experimental 
and computational studies on Al n C~ clusters. Experi- 
mental mass spectrum on Al~ show the standard magic 
number at Alj~ 3 and Al^" 3 due to electron shell closing. 
Among those anionic A1„C~ clusters, AI7C - is found to 
be particularly stable. The ground state structures of the 
A1„C _ (n — 1 — 13) clusters are determined from density 
functional optimizations. Theoretical calculations also 
demonstrate the high stability of A^C - , which might be 
attributed to the stable aluminum cage and large elec- 
tronic gap. We suggest that the magic AI7C - cluster 
could be building blocks of future cluster-assembled ma- 
terials. 
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